Chiral order of Spin-1/2 frustrated quantum spin chains 
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The ordering of the frustrated 5'= 1/2 XY spin chain with the competing nearest- and next- 
nearest-neighbor antiferromagnetic coupUngs, Ji and J2, is studied by using the density-matrix 
renormaUzation-group method. It is found that besides the well-known spin-fiuid and dimer phases 
the chain exhibits gapless 'chiral' phase characterized by the spontaneous breaking of parity, in 
which the long-range order parameter is a chirality, K;=S'f Sj^^j — Sf Sf^i, whereas the spin corre- 
lation decays algebraically. The dimer phase is realized for 0.33 <j = J2/J1 ^ 1-26 while the chiral 
phase is realized for j > 1.26. 



Considerable attention has been paid to the ordering 
of frustrated quantum spin chains since these systems 
exhibit a wide variety of magnetic phases due to the in- 
terplay between quantum effect and frustration. We con- 
sider here an anisotropic frustrated quantum spin chain 
described by the XXZ Hamiltonian, 
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where is the spin-S* spin operator at the £th site, 
Jp > is the antiferromagnetic interaction between the 
nearest-neighbor (p = 1) and the next-nearest-neighbor 
{p = 2) spin pairs, and A (0 < A < 1) represents an ex- 
change anisotropy. Note that A = and A = 1 represent 
the XY and Heisenberg chains, respectively. 

In the classical limit, 5 — > 00, the system in the ground 
state exhibits a magnetic long-range order (LRO) char- 
acterized by a wavenumber q. The order parameter is 
defined by 
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where L is the total number of spins. While the LRO is 
of the Neel-type {q — tt) when the frustration is smaller 
than a critical value, i.e., j = J2/J1 < 1/4, it becomes 
of helical-type for j > 1/4 with q = cos~^ (— l/4j). 
It should be noticed that both the time-reversal and 
the parity symmetries are broken in this helical ordered 
phase. In the XF-like case (0 < A < 1), the helical or- 
dered state possesses a two-fold discrete chiral degener- 
acy characterized by the right- and left-handed chirality, 
in addition to a continuous degeneracy associated with 
the original U{1) symmetry q£ the XY spin. The chiral 
order parameter is defined bytl 
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The chiral order parameter changes its sign under the 
parity operation but is invariant under the time-reversal 
operation. 

In contrast to the classical case S* = cx), in the quan- 
tum case 5 < 00 the magnetic LRO (||) isL-completely de- 
stroyed by strong quantum fluctuations.a Nevertheless, 
for the case of the S = 1 XY-like chain, we have recently 
shown on the basis of the exact-diagonalization (ED) 
and the density-matrix renormalization-group (DMRG) 
methods that there exist novel chiral phases in which 
only the chirality (ra) exhibits a LRO without the mag- 
netic helical LRO ^).u Two distinct chiral phases have 
been found, one gapless and the other gapped. With in- 
creasing J, the system exhibits two phase transitions first 
from the Haldane phaseQ to the gapped chiral phase, and 
then from the gapped chiral phase to the gapless chiral 
phase. These chiral ordered phases break only the par- 
ity symmetry spontaneously, with preserving both the 
time-reversal and the translational symmetries. We also 
constructed the ground, state phase diagram of the model 
(|I|) for the 5 = 1 caseMel 

By contrast, the situation is less clear in the case of 
the 5 = 1/2 XY chain. There remains a discrepancy 
between, the field-theoretical resultaJ and our numerical 
resultsu While the field-theoretical analysis based on the 
bosonization method predicted the occurrence of the chi- 
ral phase in the large j region, the Binder parameter of 
the chirality calculated via the ED method up to L = 20 
exhibited no sign of the chiral phase. Although Aligia 
et alB recently pointed out that the system size L = 20 
might be insufficient to deal with the chirality in the large 
j region, the question whether the chiral phase is realized 
in the 5=1/2 XY chain has not been clarified so far. 

In the present paper, we try to resolve this discrepancy 
by determining the phase diagram of the 5=1/2 XY 
chain (|l]) based on the DMRG method. The inethod 
used is the same as that in the previous work.Q Using 
the DMRG method, we calculate appropriate correlation 
functions associated with each order parameter charac- 
terizing the phases, and analyze their long-distance be- 
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haviors. The results of our analysis strongly suggest the 
appearance of the chiral ordered phase for j > jc(— 1-26). 

Before going into details of our calculation, we briefly 
summarize the known properties of the ground-state 
phases of the S" = 1/2 spin chain (1). The system under- 
goes a phase transition from the spin-flja|ii-phase to the 
dimer phase at j — jd with increasing j.Erll3 The critical 
value-a'd has been estimated to be jd ~ 0.33 for the XY 
case.t3 The spin-fluid phase realized at j < jd is charac- 
terized by gapless excitations above the ground state and 
an algebraic decay of spin correlations, while the dimer 
phase realized at j > jd is characterized by a finite en- 
ergy gap above the doubly degenerate ground states and 
an exponential decay of spin correlations. In the dimer 
phase, both the parity and the translational symmetries 
are broken spontaneously. The order parameter charac- 
terizing the dimer phase is given as 
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The dimer phase is further divided into two regions by the 
so-called Lifshitz line at j = j^ according to the nature of 
the spin correlations: Whereas the structure factor S{q) 
has a maximum at q = tt for j < jL, the wavenumber 
q of the maximum of S{q) shifts to an incommensurate 
value q = Q < TT for j > ji^. This incommensurate 
character is regarded-.as the vestige of the helical order 
in the classical case.E3 Meanwhile, the question whether 
the chiral phase exists for larger j for the XY chain still 
remains controversial as mentioned above. 

In order to probe the possible phases, we calculate the 
chiral, dimer and spin correlation functions defined by 
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which correspond to the order parameter (H), and 
(^, respectively. Here tq represents the center position 
of open chain, i.e., rp = L/2 for even r and tq = (iy + l)/2 
for odd r. We employ the infinite-system DMRG algo- 
rithm introduced by White.tfl The notation (• ■ •) repre- 
sents the expectation value in the lowest energy state in 
the subspace of S^^^^^ = 0. Note that, since jtiie chains 
treated in our calculation are sufficiently large|i3 one can 
safely avoid the finite-size effect arising from the incom- 
mensurate character of the spin correlation as pointed 
out by Aligia et ala 

The calculated r-dependence of the chiral, dimer, and 
spin correlation functions are shown in Fig|l| (a)-(c) on 
log-log plots for several typical values of j. As can be seen 
from Fig.|l| (a), the data of Cf^ir) for j > j^ — 1.26 are 
bent upward at larger r suggesting a finite chiral LRO, 
while for j < jc they are bent downward suggesting an 
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(b) dimer correlation 
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FIG. 1. Correlation functions versus r on log-log plots 
for various j: (a) chiral correlation dir); (b) dimer corre- 
lation ( — l)''C^ini('')i (c) spin correlation Cf{r) divided by 
the oscillating factor cos(Qr). The number of block states 
kept in the DMRG calculation is m = 450. To illustrate 
the m-dependence, we also indicate by crosses the data for 
m — 400 and 350 for several cases where the m-dependence 
is relatively large. In other cases, the truncation errors are 
smaller than the symbols. 
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exponential decay of chiral correlations. Although the 
data around the transition point suffer from the trunca- 
tion error inherent to the DMRG calculation, we estimate 
the transition point as jc — 1.261q'q3 taking account 
of the m-dependence of the data shown in the figure. 
The appearance of the chiral gtase is consistent with the 
prediction of the field thcoryQlla 

Meanwhile, as shown in Fig.^ (b), the data of C^j„^{r) 
for j < jdc — 1.26 are bent upward for larger r suggesting 
a finite dimer LRO, whereas they are bent downward for 
j > jdc suggesting an exponential decay of dimer corre- 
lations. The dimer transition point jdc is then estimated 
to be jcd = 1.26±0.01. Figurejl] (c) exhibits the spin cor- 
relation function Cf (r) divided by the leading oscillating 
factor cos{Qr). The data of Cg{r) are bent downward 
for j < jsc — 1-26 suggesting an exponential decay of 
spin correlations, while they exhibit a linear behavior for 
j > isf suggesting a power-law decay of spin correla- 
tions.Ej This behavior of spin correlations suggests that, 
as j increases, the system exhibits a transition from the 
gapped state (j < jsc) to the gapless state (j > jsc)- We 
note that the absence of magnetic (spin) LRO haSrteen 
proven rigorously for any j and for general S < ooa 

The remaining problem is the relation among jc, jdc, 
and jsc- Two possibilities seem to be allowed from our 
data, i.e., (i) jc = jdc = jsc or (ii) jc < jdc = jsc- If the 
case (i) is realized, the system undergoes only one phase 
transition at j = jc = jdc = jsc between the dimer phase 
and the gapless chiral phase with no dimer order. If the 
case (ii) is realized, on the other hand, the system under- 
goes two successive transitions on increasing j, first at 
j = jc from the dimer phase to the "chiral dimer" phase 
where both the dimer and chiral LRO's coexist with gap- 
full excitations, and then at j — jdc = jsc from the chiral 
dimer phase to the gapless chiral phase. Although rather 
large error bars of jc and jdc prevent us from determin- 
ing which of the cases is realized, our result suggests that 
the chiral dimer phase, if it ever exists, appears only in 
a rather narrow region, less than 3.2 % of jc — 1.26, be- 
tween the dimer and gapless chiral phases. (We note that 
in the S = 1 XY chain the gapped chiral phase exists for 
0.473 <j < 0.49, whose width corresponds to about 3.6 
% of jc — 0.473oQ) Further work will be necessary to 
solve the problem whether the chiral dimer phase exists. 

In summary, from numerical studies of the ground- 
state properties of the frustrated S' = 1/2 XY chain (|^), 
we have found that 

1. The chiral LRO is realized for a range of j > jd— 
1.26). This is consistent with the prediction by 
Nersesyan et. alB 

2. The dimer phase exists in the range jd{— 0.33) < 
J < jd,(~1.26). 

Detailed results including the critical properties and the 
full phase diagram in the j-A plane will be reported else- 
where. 
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In the infinite-system DMRG algorithm, the system size 
increases by two at each DMRG step and the calculation is 
continued until the L-convergence of the data is attained. 
We have performed the calculation more than 750 DMRG 
steps (corresponding to the system with L = 1500 sites) 
and confirmed the L-convergence. 

We consider that the failure of the ED analysis of Ref. ^ 
can be ascribed to the finite-size effect discussed in Ref. 
1^. The wavenumber Q characterizing the incommensura- 
bility decreases rapidly from Q — n at j — Jl to Q = n /2 
at J — > oo as j increases. In the 5* = 1/2 XY chain, the 
shift of Q from n/2 turns out to be smaller than 0.037r 
for j >jc — 1.26. Meanwhile, in a finite open chain with L 
sites, because of the condition that the wavefunction should 
vanish on both ends of the chain, the incommensurability 
smaller than 27r/L cannot be taken into account. As a con- 
sequence, the system size must be larger than L ~ 70 if 
one wishes to detect the shift of Q of order 0.037r. Thus, 
the calculation on chains up to L = 20 might fail to extract 
the true asymptotic properties for j >jc- 
The data of CJ(r) shown in Fig.|^ (c) are largely scattered, 
which might be attributed to the possible influence of cor- 
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rection terms characterized by wavenumbers Q' 7^ Q. Nev- 
ertheless, we consider that we can judge the behavior of 
the leading term of Cf (r) from the figure and make the 
conclusion in the text about (r). 
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